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ABSTRACT 

Located at z = 0.203, A2163 is a rich galaxy cluster with an intra-cluster medium (ICM) that exhibits extraordinary properties, 
including an exceptionally high X-ray luminosity, average temperature, and a powerful and extended radio halo. The irregular and 
complex morphology of its gas and galaxy structure suggests that this cluster has recently undergone major merger events that involve 
two or more cluster components. In this paper, we study the gas structure and dynamics by means of spectral-imaging analysis of 
X-ray data obtained from XMM-Newton and Chandra observations. From the evidence of a cold front, we infer the westward motion 
of a cool core across the E-W elongated atmosphere of the main cluster A2163-A. Located close to a galaxy over-density, this gas 
'bullet' appears to have been spatially separated from its galaxy (and presumably dark matter component) as a result of high-velocity 
accretion. From gas brightness and temperature profile analysis performed in two opposite regions of the main cluster, we show 
that the ICM has been adiabatically compressed behind the crossing 'bullet' possibly because of shock heating, leading to a strong 
departure of the ICM from hydrostatic equilibrium in this region. Assuming that the mass estimated from the Y x proxy best indicates 
the overall mass of the system and that the western cluster sector is in approximate hydrostatic equilibrium before subcluster accretion, 
we infer a merger scenario between two subunits of mass ratio 1:4, leading to a present total system mass of M 500 = 1.9 x 10 15 Mo. 
Additional analysis of the spatially-separated northern subcluster A2163-B does not show any evidence of strong interaction with the 
main cluster A2163-A, leading us to infer that the physical distance separating the northern subcluster and the main component is 
longer than the projected separation of these components. The exceptional properties of A2163 present various similarities with those 
of 1E0657-56, the so-called 'bullet-cluster'. These similarities are likely to be related to a comparable merger scenario. 
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1. Introduction ture has been seen in weak lensing (see e.g. lSquires et al.ll 1997b 

iRadovich et alj 20081) and in the galaxy density distribution 

A2163 is a rich galaxy cluster (richness class 2) located at (Maurogo rdato et al.ll2008l). Moreover. A21 63 is known to host 

z = 0-203. AfteMnitial X-ray detection by HEAO 1 A-l a prominent radio halo (iFeretti etZll200ll 12001) and to be a 

dKowalski et alj [1984J), combined observations by the Ginga possible s ource of non t hermal hard X-ray emission, as sug- 

mH , and Einstein X-ray satellites revealed the extraordinary prop- ge sted bv iRephaeli et all dUS) from RXTE observations and 

^ • ernes of its hot gas content, with exceptionally high lumi- Mmion & Alien! (l2009h from Chandra data analysis. Due to its 

H , nosity and average temperature (L x = 3.5 x 10 JS W, kT = high mermal energy content , A2163 is also a favoured targe t 

, 13.9 keV, | Arnaud et alj 112230, suggesting a very high clus- for sz obser vations. First detected by Wilba nks et all (Tl 994) 

ter mass. Extensive follow-up observations have further re- and used t0 det e rmi ne the Hubble constant in iHolzanfel et al l 

vealed several signatures of major merger events at van- (fl997l), A2163 is the first galaxy cluster detected in the submm 

ous wavelengths in this hottest Abell cluster. ROSAT obser- band jLamarre et al.lfl998l A recent analysis of its SZ signal bv 

vations showed an irregular X-ray morphology dElbazetalJ | N ord et al.1 d200lT nas confirmed the exceptionally high temper- 

[19951), while ASCA observations revealed complex thermal ature of its hot gas content independently of X-ray analyses, 
structure (Markevitch et al. 1994) , more recently confirmed 

from Chandra data dMarkevitch & Vikhlininl200lt iGovoni et alj Quite apart from the evidence for a very high total system 

l2004t lowers et al] f2009). Evidence of a clear cluster substruc- mass, the exceptionally large amount of thermal energy trapped 

within the hot atmosphere of A2163 is likely to be related to on- 

Send offprint requests to: H. Bourdin going cluster merger events. For this reason, the mass of A2163 

1 10 38 W = 10 45 erg s 1 ; X-ray luminosity in the 2-10 keV band has been challenging to estimate from X-ray analyses, which 

has been corrected for luminosity distance assuming H = 70 km s 1 rely on the gas being in hydrostatic equilibrium. Indeed, de- 

Mpc -1 , and A = 0.7. parture of the cluster gas content from hydrostatic equilibrium 
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Fig. 1. XMM-Newton -EPIC exposure mosaic (left) and Chandra-ACIS exposure (right). 



may explain t he 50 % dispersio n obs erved between mass esti - 
mates of e.g. Elba z et aQ (1 19951) and iMarkevitch et ail (|l996), 
performed assuming different temperature gradients in the clus- 
ter outskirts. It is also interesting to note the discrepancy between 
X-ray measurements of the gas tempe rature and the i ndirec t 
weak lensing temperature estimate of ICvpriano et a l. (2004), 
kTsis = 6.63^°g keV, obtained assuming energy equipartition 
between galaxies and gas. 

Our recent analysis of opt i cal dat a (WFI/2.2m, VIMOS/VLT, 
UT3) in Maurogor dato et all (120081 hereafter Paper I) revealed 
the complex galaxy structure and dynamics in A2163, where 
a main cluster, A2163-A, is connected to a smaller subcluster 
to the North, A2163-B, via a galaxy bridge. The main cluster, 
A2163-A, itself has two brightest galaxies (BCG1 and BCG2), a 
bimodal morphology, and an exceptionally large velocity disper- 
sion (1434 ± 60 km s~'). As part of a multi-wavelength analysis 
programme aimed at investigating the impact of cluster mergers 
on galaxy and ICM evolution, the present article focuses on the 
analysis of XMM-Newton and Chandra observations of the ICM 
in A2163. This work particularly aims to investigate the merger 
scenario involving the various components of A2163, and to 
provide new insights on the merger impact on X-ray mass esti- 
mates. Data preparation and analysis issues related to the XMM- 
Newton and Chandra observations are detailed in Sections |2]and 
[3] Results obtained on ICM substructure and merger dynamics 
are discussed in Section [4] while Sections and [6] concern the 
gas pressure structure and mass of A2163. Discussion and con- 
cluding remarks are provided in Section [7] Except if otherwise 
noted, the confidence ranges on individual parameter estimates 
are 68 %. In the following, Mg is defined as the mass within the 
radius ig at which the mean mass density is 6 times the critical 
density of the universe at the cluster redshift. Moreover, intra- 
cluster distances are computed as angular diameter distances, 
assuming a A-CDM cosmology with Ho = 70 km s _1 Mpc -1 , 
CIm = 0.3, Q.a = 0.7. Given these assumptions, an angular sep- 
aration of 30 arcsec corresponds to a projected intra-cluster dis- 
tance of 100 kpc. 



2. Observations and data preparation 

2.1. Observations 

2.1 .1 . The XMM-Newton data set 

The XMM-Newton data set is a mosaic composed of five ob- 
servations obtained with the European Photon Imaging Camera 
(EPIC), and available in the XMM-Newton archive (observa- 
tion IDs: 01 12230601, 01 12230701, 01 12230801, 01 12230901, 
01 12231001). This data set has been used to map the ICM ther- 
mal structure in the central region of A2163, and to investigate 
its radial brightness and temperat ure structure on a large scale . 
Following a similar scheme as in lBourdin & Mazzotta (2008), 
hereafter BM08, we filtered the data set through spatial and tem- 
poral wavelet analyses in order to remove the contribution of 
point sources and soft proton flares. A summary of the effective 
exposure time remaining after filtering is provided in Table Q] 

2.1 .2. The Chandra data set 

The Chandra data set is a 72 ks pointing observation of A2163 
obtained with the AXAF-I CCD Imaging Spectrometer (ACIS- 
I) and available in the Chandra archive (observation ID: 01653). 
These data have been used to map the gas brightness structure at 
high angular resolution near the cluster centre. Following a sim- 
ilar procedure as for the XMM-Newton data, we filter events in 
order to remove point source and soft proton flare contributions 
(see effective exposure in Tab.[T]i. 

2.2. Data preparation 

In order to perform self-consistent brightness measurements and 
spectroscopic estimates using the XMM-Newton or Chandra 
data sets (see Section |3J, we re-sampled photon events within 
fixed grids matching the angular and spectral resolution of 
each focal instrument in sky coordinates (k, I) and energy (e). 
Events associated with the five XMM-Newton pointings have 
thus been gathered into three event cubes corresponding to each 
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Table 1. Effective exposure time of each XMM-Newton -EPIC and Chandra-ACIS observation. In brackets: fraction of the useful 
exposure time after solar-flare "cleaning". 



XMM-Newton 


Centre coordinates 


MOS1 effective 


MOS2 effective 


PN effective 


obs. IDs 




exposure time (ks) 


exposure time (ks) 


exposure time (ks) 


0112230601 


16hl5m46.01s -06°09'00.0" 


8.8 (56.2 %) 


9.8 (63.5 %) 


5.7 (56.2 %) 


0112230701 


16hl6m48.00s -06°09'00.0" 


11.8(42.5 %) 


12.1 (43.4%) 


6.1 (29.9%) 


0112230801 


16hl5m46.01s -05°51'00.0" 


16.1 (57.5 %) 


14.6 (52.2%) 


10.2 (49.0 %) 


0112230901 


16hl4m34.01s -06°09'00.0" 


7.4 (27.5 %) 


6.9 (25.5 %) 


3.3 (16.7 %) 


01122301001 


16hl5m46.01s -06°27'00.0" 


19.5 (69.3 %) 


16.5 (58.8 %) 


15.4 (72.4 %) 


Chandra obs ID 


Centre coordinates 


ACIS effective exposure time (ks) 


01653 


16hl5m45.77s -06° 08'55.0" 




65.2 (91.7%) 





of the three EPIC instruments, while events associated with 
the Chandra pointing have been binned into a unique event 
cube. Following a similar procedure as described in BM08, we 
then associated a set of local 'effective exposure' and 'back- 
ground noise' arrays to these event cubes. The effective exposure 
E(k, I, e) is computed as a linear combination of exposure times 
associated with each pointing, with correction for spatially vari- 
able mirror effective areas, filter and other focal instrument trans- 
missions, CCD pixel area with correction for telescope motion, 
gaps and bad pixels^. The 'background noise' array (k, I, e) has 
been modelled as the sum of 3D functions accounting for galaxy 
foreground and Cosmic X-ray Background (CXB) emissivities, 
(<?), and also false photon detection from the particle induced and 
out-of-time events, F p (k, I, e) and 0(k, I, e), respectively^ 

nt,ck(£, t)(k, I, e) = E(k, I, e) x n ga i ;CX B (e) 
+n ool (k, /) 0(k, I, e) 

+n p (M)F p (e). (1) 

As described in iBourdin et alj (12004), we model the parti- 
cle induced background F p (e) as the sum of a dual power-law 
continuum and a set of detector fluorescence lines (see light 
blue curves on Fig. This model has been fitted to a set 
of observations obtained during 'closed' periods of the XMM- 
Newton telescopes, and stowed periods of the Chandra ACIS 
detector. For both the XMM-Newton and Chandra observations, 
the galaxy foreground and CXB emissivities, (<?), have been es- 
timated within an external annulus located away from the target 
emission using offset pointings of the EPIC-XMM-Newton mo- 
saic. Further details concerning our modelling of the extended 
foreground and background emissions are provided in Section 

EH 

3. Data analysis 

3. 1 . X-ray brightness measurements and spectroscopy 

The event, 'effective exposure', E{k, I, e), and 'background 
noise', (k,l,e), arrays have enabled us to perform ICM bright- 
ness and average temperature measurements in various areas of 
the field of view by gathering pixels (k, 1). 

We computed brightness estimates l, x (k, I) from photon 
counts n evt (fc, /, e) in the 'soft' energy band (0.5-2.5 keV) as 
followJ3 



2 Information about these instrumental effects have been obtained 
from the XMM-Newton -EPIC Current Calibration Files (CCFs) and 
Chandra-ACIS Calibration database (CALDB 4.1) 

3 In Eq. O and {2}, spectral models are decomposed into normalised 
energy distributions, F,(e), and counts per pixel n, 
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Fig. 2. Background spectrum observable in the A2163 out- 
skirts (r>3Mpc) from the four external pointings in Table [TJ 
Light blue: particle background noise. Cyan blue: Cosmic X-ray 
Background emission. Blue and Violet: TAE emission (kT i = 
0.099 keV, kT 2 = 0.2 48 keV, see lKuntz & Snowdenll2000l and 
details in Sect. I3.3.4t . Orange: Local Hot Bubble and North 
Polar Spur residue. Red and black: overall fit and data set. 



£ x (fc,0 = 



Ze n e vt(£, e) - n hck (k, I) £ e (k, I, e) 
2 e Fi CM (kT , Z , N H , , e)E(k, I, e) 



(2) 
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where FicM(kT ,Z ,NH, ,<?) is an emission model asso- 
ciated with average source temperature, kT , metallicity, Z , 
and Galactic hydrogen absorption Nh, - Hereafter, the source 
emission spectrum Ficivi(kT, Z, Nh, e) assumes an Nh absorbed 
emission mod elled from the A strophysical Plasma Emission 
C ode (APEC. | Smith et alj|2001l) . with the element abundances 
of iGrevesse & Sauvall (1 1998b and neutral hydrogen absor ption 
cross sections o flBalucinska-Church & McCammonld 19921) . It is 
altered by a local energy response matrix M(k,l,e,e') computed 
from response matrixes files (RMF) tabulated in detector coordi- 
nates in the XMM-Newton -EPIC and Chandra- ACIS calibration 
data bases Q. 

Following this emission model, average temperatures 
kT(fc, /) are computed from spectral fitting to the emission spec- 
trum n evt (fe, l)F evt (k, I, e), registered in the energy band (0.7-12 
keV): 



T(r) = T e 



(r/r t r a 



[1 + {rlr t ) b Y' b ' 



(7) 



In Sect. 14.21 the ICM emission measure and temperature pro- 
files have been modelled by two step distributions with common 
jump radius r,, in order to calculate the physical characteristics 
of a cold front in the cluster atmosphere: 



[n p n e ] cf (r) 



D 2 n n 2 ( r /rjr 2 " 
n 2 (r/ rj r 2 " 



l+( r /r £ ,) 2 
l+(rj/r c2 )' 
1+0/ra) 2 



3/3, 



r < rs 



r > r; 



T c /(r) 



T , r<rj 
D T T ,r> rj 



(8) 



(9) 



n e vt(M)F evt (kT,Z,N H ,e) 

E(k, I, e) x n lCM (k, F ICM (kT, Z, N H , e) 

+n bc k (k, I, e). 



(3) 



3.2. X-ray brightness, ICM density and temperature profiles 

Through the above analysis, X-ray brightness and temperature 
profiles have been extracted assuming uniform emissivity within 
concentric annuli. We derived the radial brightness £ x (r) and 
associated Gaussian fluctuation cr^ir) by averaging the local 
brightnesses £ x (fc, I) of Eq. (f2j) in each profile annulus composed 
of N pixels (k, 1), as follows: 



k,l 



Y, e n e vt(^, /, <?) 



N \u [Z,F IC M(kT ,Z ,NH, , e )E(M, e )f 



(4) 
(5) 



Radial temperatures kT(r) and associated confidence inter- 
val 6kT(r) have been computed within each annulus by fitting a 
uniform emission model to the data set. To do so, we averaged 
emission models of Eq. 0) associated with each pixel (k, I) of 
the annulus, and fitted an average set of spectroscopic parame- 
ters kT(r), Z(r), Nh(V) by means of ax 2 minimisation. 

These brightness and temperature profiles have been used to 
model the underlying density and temperature of the ICM, as- 
suming spherical symmetry of the cluster atmosphere. This was 
undertaken by projecting and fitting parametric distributions of 
the radial emission measure, n p n e , and temperature, T(r) , to the 
observed profiles. Introduced by Vik hlinin et al.l d2006l) . these 
parametric distributions are: 



[n p n e ](r) = n c 



(r/r c T a 



1 



[l+(r/r e ) 2 ]3^/ 2 [l+(r/ rs ) 3 ] e/3 



'02 



[1 + (r/r c2 )2]3A ' 



(6) 



4 EPIC response matrixes are computed from canned RMFs 
corresponding to the observation period provided by the XMM- 
Newton Science Operation Centre. ACIS responses matrixes have 
been computed using the Chandra Interactive Analysis of Observations 
(CIAO) software and Chandra-ACIS Calibration database (CALDB 
4.1) 



3.3. Imaging and spectral-imaging 
3.3.1. Imaging 

The ACIS -Chandra data set has enabled us to analyse the ICM 
brightness structure in A2163 at the highest angular resolution 
available. By means of two types of wavelet analysis, we first 
mapped the ICM brightness of the large-scale cluster emission, 
then isolated the brightness structure of the central region from 
contributions related to the cluster outskirts. 

To map the large-scale ICM brightness, we first corrected 
'soft' photon (0.5-2.5 keV) counts for vignetting and back- 
ground emissivity as detailed in Eq. (0, then analysed the 
corrected image following the redund ant 'a-trous' algorithm 
dHolschneider et alJ 1989; Sh ensalll992h implemented with B3- 
spline wavelets. A de-noised brightness map was then con- 
structed from the overall wavelet transform of the image, 
with coefficients thresholded iteratively assuming that photon 
counts follow a Poi sson distribution (see [giiaoui & Rue 1995; 
IStarck & Pierrd[T998L and references therein for details). In the 
top-left panel of Fig.|U this map is overlaid on an /?-band optical 
image obtained at the WFI/MPI 2.2m telescope (see also Paper 
I). 

To further isolate the ICM brightness structure in the clus- 
ter centre, we restricted a B3-spline wavelet analysis of a cor- 
rected photon image extracted in the [0.5-7.5] keV energy band 
to four high resolution scales, corresponding to a distance range 
of [1 - 120 arcsec]. The resulting high resolution residue of the 
ICM brightness map shown on Fig. [4] has been constructed by 
summing details of this restricted wavelet transform, with coef- 
ficients thresholded to a 4 cr confidence level. 



3.3.2. Spectral-imaging 

In order to map the ICM temperature in A2163, we used the 
EPIC-XMM-Newton data set and applied the spectral-imaging 
algorithm detailed in Bourdi n et al.l (|2004|) and B08. Following 
this algorithm, a set of temperature arrays kT(£, /, a) with asso- 
ciated fluctuations cr^ik, I, a) are first computed on various anal- 
ysis scales a, then convolved by complementary high-pass and 
low-pass analysis filters in order to derive wavelet coefficients. 
The wavelet coefficients are subsequently thresholded accord- 
ing to a given confidence level in order to restore a de-noised 
temperature map. Here, the signal analysis have been performed 
over 6 dyadic scales within an angular resolution range of 6a 
= [3.4 - 220] arcsec. This was undertaken by averaging the 
emission modelled by Eq. <[3j within overlapping meta-pixels 
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Fig. 3. Left: Galactic neutral hydrogen density column estimated in the neighbourhood of A2163, with contours of the [0.5-2] keV 
XMM-Newton -EPIC emission overlaid. Right: Correlation factor, relating IR emissivity to the Galactic Nh in a 5 square degree 
region of the sky centred on A2163. IR emissivity and Galactic Nh have been corrected for dependence on galactic latitude following 
the analytic approach of iDesert et alj d 1988b . A correlation factor p£ is computed between the average residual IR emissivity and 
Galactic Nh within sky regions corresponding to a fixed Nh interval of width £Nh = 10 21 m~ 2 . 



(k, I, a), and computing the kT(k, I, a) and cr^k, I, a) arrays by 
means of a likelihood maximisation assuming the spatially vari- 
able Galactic Nh and foreground emission model discussed be- 
low in Sect. l3.3.3l The resulting ICM temperature map shown in 
Fig. |4] was then obtained from a B2-spline wavelet analysis (see 
B08 for details) with coefficients thresholded to the 1 cr confi- 
dence level. 



3.3.3. Absorption by Galactic neutral hydrogen 

Located close to the galactic plane, in the direction of the galac- 
tic centre (1=19.1°, b=37.2°), A2163 is placed behind a rela- 
tively strong column density of galactic neutral hydrogen, Nh — 
1.5 x 10 21 crrT 2 . In this high Nh regime, temperature measure- 
ments are highly sensitive to spatial variations of the Galactic Nh 
on angular scales smaller than the cluster size, so that a precise 
modelling of these variations is required to map the ICM tem- 
pe rature. Indeed, as formerly shown by a ROSAT data analysis 
of lElbaz et al ] (1 19951) . the combination of the high average tem- 
perature of the ICM in A2163 with this high (and variable) Nh 
value makes temperature measurements highly degenerate with 
Nh estimates. For these reasons, we decided to model a priori the 
spatial variations of neutral hydrogen Galactic absorption across 
the cluster field of view. 



molecular cloud 0, the global budget of Hydrogen density 
column is dominated by neutral hydrogen in this region of 
the sky, enabling us to use galactic dust emissivity maps 
at 100 //in as a tracer of Nh spatial variations. To model 
these variations, we follow a two step approach where large 
scale variations related to changes in the Galactic Interstellar 
Radiation Field are separ ated from smaller scale fluctuations . 
We used the IRAS/IRIS dMiville-Deschenes & Lagachdf2 005) 
and Leiden/ Argentine/Bonn (LAB. iKalberla et alJ 120051) sur- 
veys of the all sky IR emissivity and Galactic Nh, and corrected 
their dependence with galactic latitude, b, by subtracting two 
cosecant laws with shap e -^n , following the analytic approach 
of iBoulaneer & Peraulj (119881) and equations (3.1) and (3.2) of 
IDesert et al.1 0988). We then measured the large-scale correla- 
tion factor, relating the residual IR emissivity and Galactic 
Nh in a 5 square degree region of the sky centred on A2163 (see 
also Fig. [3]). We finally model the smaller scale spatial variations 
of Nh across the cluster field of view, by applying the ^ factor 
to the IRAS/IRIS maps, and deducing the expected Nh map at 
high angular resolution. As shown in Fig. [3] this map reveals us 
Nh variations of the order of 20 %, with a NorthEast-SouthWest 
gradient that could have strongly distorted X-ray spectroscopic 
estimates, if not taken into account. 



3.3.4. Foreground and background spectral contributions 



To investigate the spatial variations of Galactic Hydrogen 
density column up to an angular resolution of about one ar- 
cminute, we used Galactic dust emissivity maps at 100 fim. 
The spatial variations of the overall Hydrogen (i.e., both neu- 
tral HI, and molecular H2) density column are indeed spa- 
tially correlated with g alactic dust emissivity, ttr, in the 100 
I R band (see e.g. iBoulanger et alj [1996: Snowden et alj 
119981) . Moreover, as A2163 is located outside any Galactic iDame et aiTd2001h . 



The spectral 'background noise' model nt, c k (k, I, e) includes not 
only contributions related to the instrumental noise, but also 
those due to the extended emission of Galactic foregrounds 
and the Cosmic X-ray background. Following the approach de- 
scribed in BM08, we model the Cosmic X-ray background with 



See e.g. the large-scale CO survey of the Galactic plane by 
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Fig. 4. A2163. Left panels: X-ray e mission overlaid on galaxy maps. Optical observations performed at the MPI 2.2m telescope, see 
Paper I dMaurogordato et aUl2008l) . X-ray images obtained from wavelet analyses of a Chandra exposure in the 0.5-2.0 keV band. 
Top-left: overall cluster emission. Bottom-left: high resolution analysis of the ICM emission near the cluster centre, see Section [331 
for details. Right panels: ICM temperature maps obtained from wavelet spectral-imaging analysis of the EPIC-XMM-Newton data 
set. Top-right: black iso-contours from the Chandra image (same as top-left map) overlaid as black iso-contours. Bottom-right: black 
iso-contours from high resolution Chandra residue (same as bottom-left map); white iso-density contours from galaxy density map 
(see also Fig. 6 in Paper I). 



a power law of index y — 1 .42 (see e.g. iLumb et aDl2002l) and 
the Galactic foregrounds by the sum of several absorbed and un- 
absolved thermal components. 

As already not ed by iPratt et alJ (l200lh and 

Markevit ch & Vikhlininl d2001l) . the Galactic foreground 
emission near A2163 includes a strong excess in the soft band 
(E < IkeV), with regards to the average emission expected 
from 'blank sky' observations. Taking into account this excess, 
in addition to the local variations of Nh (see Section 13.3.3b 
is critical to perform any spectral analysis, in particular when 
measuring th e ICM temperatu re in the cluster outskirts. As 
suggested by IPratt et all (l200lh . the origin of this soft excess 



is probably related to the location of A2163 near the North 
Polar Spur (NPS). For this reason, we decided to model 
the Galactic foreground emission by the sum of an unab- 
solved thermal emission accounting for the local hot bubble 
(LHB, YYlhb = 0.1 keV, see e.g. iKuntz & Snowd en 2000) 
and three absorbed thermal components accounting for the 
Galactic 'transabsorptio n emission' (TAE ; kTi = 0.099 keV, 
kT 2 = 0.248 keV, see IKuntz & Snowderj|2000h and the NPS 
(see Fig. [2j. In order to constrain this complex spectral model, 
we took advantage of the sky coverage provided by the four 
offset pointings available in the XMM-Newton archive (see 
Table [TJ and fitted the various background components to the 
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Fig. 5. Top profiles: ICM brightness (left) and projected temperature (right) profiles corresponding to the sector illustrated in the 
bottom-left panel of Fig. |4] Bottom profiles: ICM density and 3D temperature modelled as disrupted distributions with common 
jump position (the dot-dashed line, see also Sect. l4.2l and Eq.[6]and[7]l. The ICM emissivity corresponding to these radial distribution 
has been fitted to the ICM brightness and projected temperature, as shown by dotted line on the top profiles. 



X-ray emission from an external annulus (2700 < r < 4000 
kpc) where no cluster emission is expected to be detectable- 
Modelling the spatial variations of Galactic Nh as detailed 
in Section I3.3.3I and fixing the temperature of the thermal 
components representing the galactic TAE and LHB enabled 
us to constrain the emissivity and temperature of the excess 
component. The temperature value obtained (kTvps - 0.3keV) 
is in agreement wit h expectations for emi ssion from the North 
Polar Spur (see e.g.lwillingal e et al.ll2003l) . 



4. ICM dynamics and merger events 

4.1. Gas distribution and thermal structure of A2163 

The Chandra-ACIS brightness map and iso-contours in the top 
left and right panels of Fig. [4] show us the complex and irreg- 
ular morphology of the ICM in A2163. The map and contours 
reveal a main component (A2163-A) and a Northern subcluster 
(A2163-B), the main component itself appearing irregular with 
an eastern extension on a large scale and a triangular shape in 
the central region. The high resolution residue obtained from 
wavelet analysis of the Chandra-ACIS 'soft' image -as detailed 
in Sect. I3.3l and shown in the bottom-left panel of Fig. |4]- en- 
ables us to resolve further details in the central region, and to 
separate a brightness peak located close to the brightest cluster 
galaxy (BCG1, see also Paper I) from a sharp-edged, wedge- 
shaped feature located to its west. 



6 This radial range corresponds to a region located well beyond r 2 oo, 
as estimated from optical measurements in Paper I 



The XMM-Newton -EPIC temperature map in Fig.|4]shows 
the very hot mean temperature and anisotropic thermal struc- 
ture of the ICM near the centre of the main cluster A2163-A, 
and also the colder temperature (kT ^ 6keV) of the Northern 
subcluster, A2163-B. The map exhibits in particular a hot re- 
gion to the North of A2163-A (kT a 18keV) and a colder 
clump to the SouthWest (kT =s lOkeV). This overall thermal 
struct ure is consi s tent with former analyses of C ha ndra data by 
e.g.lGovoni et al.ld2004l) . lMilrion & Allen! d2009l) or lOwers et alj 
(I2009I) . apart from slight morphological differences in the shape 
and location of the hot region in our maps. This hot region is 
indeed spatially coincident with a prominence of Galactic IR 
emissivity, likely leading to a region of enhanced absorption (see 
Fig. |3j. The better sensitivity of the XMM-Newton mosaic fur- 
ther enables us to detect a clear drop of the ICM temperature 
toward the outskirts of A2163-A (r > 700 kpc). 

4.2. A gas 'bullet' crossing the main cluster atmosphere 
separately from galaxies 

In the bottom-right panel of Fig. |4] we show the high resolu- 
tion residue obtained from the Chandra analysis overlaid on the 
XMM-Newton -EPIC temperature map. Interestingly, we notice 
that the Southwestern cold clump in A2163-A is associated with 
the wedge-like feature revealed from the Chandra data. The cold 
clump thus appears as a cool core embedded in the hotter ICM 
of A2163-A, that is delimited in the westward direction by a 
brightness edge. 

In order to investigate 3D structure of this cool core, we ex- 
tracted the brightness and temperature profiles corresponding to 
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an ICM sector incorporating the brightness edge, as illustrated 
in the bottom-left panel of Fig. [4] As shown on Fig. [5] these 
two profiles have enabled us to model the ICM density p(r), 
and 3D temperature T(r), as two disrupted functions defined in 
Eq. (O and (O, with common jump location and amplitudes. 
The density and temperature jump amplitudes thus measured are 
D n = 21 = 1-28!^ and D T = £ = 1.65*$, respectively, 
where indices 1 and 2 refers to internal and external core re- 
gions. Consistent with continuous pressure at the jump location 
(Dp = D„Dt = 1.29+0 ^3), these jumps are likely to be associ- 
ated with a cold f ront. This co l d fron t has also been seen in the 
same data set by lOwers et al.l d2009l) . who measured a similar 
pressure jump to ours from analysis of a region corresponding to 
the Southern part of our sector where the density jump appears 
as more prominent in project ion. By analogy with 1 E0657-56 
(the so-called 'bullet-cluster', M arkevitch et al.l l2002). the loca- 
tion of the cold front with respect to the wedge-shaped residual 
emission might indicate the westwards motion of a stripped cool 
core embedded in the hotter atmosphere of A2163-A. 

To try and constrain the history of a possible subcluster in- 
fall associated with this cool core, in the bottom-right panel of 
Fig. [4] we overlay the galaxy iso-density contours obtained from 
WFI data (Paper I) on the ICM temperature map of A2163. 
As discussed in Paper I, the galaxy iso-density contours map 
a complex, bimodal galaxy distribution in A2163-A. This bi- 
modal distribution corresponds to an E-W elongation of the clus- 
ter dark matter halo, a s revealed from weak lensing analysis of 
Rado vich et al.l d2008l) . Interestingly, we now observe that the 
less dense of the two galaxy subclusters in A2163-A is located 
near the cool core, but is separated from it at a projected dis- 
tance of about 30 arcsec. Its location to the west of the cool core 
suggests a scenario in which a subcluster has crossed the A2163- 
A complex from east to west; the offset wedge-shaped residual 
emission then suggests subcluster gas loss gas due to ram pres- 
sure of the main cluster atmosphere, leading to its present sepa- 
ration from the constituent galaxies. 



4.3. The subcluster A2163-B 

4.3.1. A2163-B. Gas properties and mass 

Located at close projected distance (=; 1.25 Mpc) from the main 
cluster A2163-A, A2163-B appears in Fig.|4]as a small subclus- 
ter with regular morphology, hosting colder gas than A2163-A. 
Due to the small angular distance separating A2163-B from its 
bright companion, the X-ray emission from both A2163-A and 
A2163-B must be analysed simultaneously. We separated the 
two cluster emissions by modelling the emissivity from A2163- 
A in each energy band from profiles p(r) and T(r), used to model 
its gas density and temperature structure in Sect. [6] To do so, 
a 2D yS-model was first fitted to the surface brightness distribu- 
tion of A2163-A, enabling us to derive ellipticity parameters and 
to correct the projected emissivity, I.(d,e) = 2 J^°° rdr, 
for these parameters. In order to further reduce the contribution 
from A2163-A, the spectra from A2163-B have been extracted in 
a northern semi-sector centred on the subcluster emission peak 
and extending away from the main cluster emission. 

This spectral component separation has enabled us to de- 
rive the residual brightness and temperature profiles of A2163- 
B. As shown in Fig. [6] the gas temperatures obtained are lower 
than expected from visual inspection of the temperature map of 
Fig- 111 since here the emission contribution from the main clus- 
ter has been removed. Fitting this data set with analytical den- 
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Fig. 6. A2163-B. Top: ICM surface brightness profile. Light 
blue: blank-sky background emissivity. Blue: blank-sky back- 
ground and main cluster emissivity. Middle: projected ICM tem- 
perature profile. Bottom: Mass profile modelled from the surface 
brightness and projected temperature profiles, assuming ICM 
hydrostatic equilibrium. The ICM surface brightness and pro- 
jected temperature have been modelled analytically (see dotted 
lines on the corresponding plots). 



sity and temperature profiles has enabled us to derive a mass 
profile of A2163-B, assuming hydrostatic equilibrium. The clus- 
ter gas and total mass estimates derived from this analysis are 
M 



500 = 0.22+° : °2 x 1O 14 M and M500 — xu.^™ ai@t 
spectively, where confidence intervals have been obtained from 
random realisations of parameters in the fitted density and tem- 
perature profile of Eq. © and (0. 



2.1 + 0.1xl0 14 M o ,re- 



4.3.2. A2163-B. interaction with main cluster, A2163-A 

Using galaxy velocity measurements from paper I and applying 
the two-body formalism, we calculated a set of bounded solu- 
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Fig. 7. A2163 gas temperature map with two complementary 
profile extraction regions overlaid. The blue region is expected 
to be located ahead of the moving cool core, while the red region 
is expected to be located behind. Both regions are centred on the 
X-ray emission peak. 

tions for the dynamics of the merger system A2163-A - A2163- 
B, with various separation distances including close solutions 
where both subclusters are almost observed in the sky plane, 
and more distant solutions where both clusters are observed as 
aligned along the line of sight (see Appendix for details). The 
ICM properties of A2163-B do not show any strong evidence of 
interaction with the main cluster atmosphere. In particular, the 
gas brightness and temperature maps in Fig. [4] do not exhibit 
any strong thermal or pressure enhancements in the direction of 
A2163-A within an over-density radius of rsoo, in contrast with 
expectations of hydrody namical N-body simulat ions of bimodal 
cluster mergers (see e.g. Ricker & Sarazin 2001). Moreover, the 
gas mass fraction in A2163-B is close to the value expected for 
an isolated system (f gas ,5oo — 0.1). From the lack of any close 
interaction evidence with A2163-A, it is likely that A2163-B is 
physically separated from its companion at a distance larger than 
the radius, rsoo, of the main cluster. 

5. Merger effect on the gas pressure 

The evidence for a gas 'bullet' separated from its galaxies cross- 
ing the A2163 atmosphere suggests that the main cluster has re- 
cently accreted a subcluster along the East-West direction (see 
Section I4.2l i. In order to investigate the effects of this major 
merger event on the ICM thermodynamics, we extracted three 
sets of ICM brightness, temperature and pressure profiles cor- 
responding to the eastern and western cluster sectors shown on 
Fig. [7] and to the overall cluster A2163-A. 

5.1. Brightness and temperature profiles 

Extracted from the XMM-Newton data set and shown on Fig. [8] 
the ICM brightness and projected temperature profiles have been 
centred on the eastern X-ray peak revealed on the Chandra high 
resolution map of Fig. [4] and located close to the cluster BCG1 
as shown in Paper I. This choice has enabled us to centre the 
profiles near the expected cluster mass centre, while radially 



integrating the innermost gas brightness distribution along the 
X-ray isophotes. The maximum radius of the profiles has been 
roughly set to a distance where cluster reaches the 'background 
noise' emissivity in the soft X-ray band. Due to the strong East- 
West asymmetry of the cluster atmosphere we are able to extract 
profiles up to larger radius on the eastern side. 

We first observe that the shapes of the eastern and west- 
ern brightness profiles differ strongly near the centre of A2163- 
A, due to the asymmetry introduced by the crossing 'bullet'. 
The two profiles cross each other at a radial distance of about 
3.5 arcmin, beyond which the eastern brightness profile exceeds 
the western one. As expected, the overall cluster profile aver- 
ages the radial brightness values corresponding to the eastern 
and western cluster sectors. Given the limited statistics avail- 
able, the eastern and western temperature profiles show only 
marginal differences, most noticeably for the outermost radial 
bins and in the second radial bin intercepting the western cool 
core. Interestingly, we observe that the overall temperature pro- 
file seems to reflect the higher temperatures on the eastern cluster 
side, where the radial ICM emissivity is the highest. 



5.2. Pressure profiles 

As shown in the bottom-left panel of Fig. [8] the ICM brightness 
and temperature profiles have been used to model three pres- 
sure profiles corresponding to the eastern, western and overall 
cluster outskirts (r > 400 kpc). The ICM density and 3D tem- 
perature profile have been modelled following parametric radial 
distribution functions -see Eq. © and and multiplied by 
each other to yield the pressure. As observed for the brightness 
profiles, the eastern and western pressure profiles differ clearly 
in the innermost region of A2163-A, then cross each other at a 
radial distance of about 800 kpc, above which the eastern profile 
has a higher pressure than the western profile by a factor ~ 2. 
This pressure excess may be related to ICM shock heating be- 
hind the crossing 'bullet', leading to an adiabatic compression 
of the gas. The overall cluster pressure profile seems to reflect 
the higher pressure values corresponding to the eastern cluster 
sector, so that the ongoing merger event is likely to have strongly 
disturbed the shape of the overall cluster pressure profile. 

A recent analysis o f a representative sam ple of nearby galaxy 
clusters (REXCESS, iBohringer etal]|2007h has shown that the 
ICM pressure is expected to follow an average profile with low 
dispe rsion partially relat ed to the thermodynamical state of clus- 
ters dArnaud et al.ll20l61 hereafter A10). In order to be compared 
with average cluster pressure profile of the REXCESS sample, 
the ICM pressure profiles of A2163-A have been rescaled to 
characteristic pressures P500 related to the cluster mass M500, 
in the framework of self-similar cluster evolution models (see 
Eqn. 5 of A10). As detailed in the following Section, the clus- 
ter mass we used for this scaling was estimated from the Yx 
proxy. The bottom-right panel of Fig. [8] shows that the scaled 
profiles associated with the overall cluster and the western sector 
both follow the universal shape of the REXCESS cluster sam- 
ple, while the scaled profile associated with the eastern sector 
appears shallower. Consistent with the discrepancy observed by 
A10 between relaxed and morphologically disturbed clusters in 
the REXCESS sample, this trend favours a scenario where the 
western side of A2163-A shows the pressure structure of a self- 
similar cluster, while the pressure structure in the eastern cluster 
side has been strongly disturbed by the ongoing merger event. 
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Fig. 8. A2163-A profiles associated with the overall cluster (in black) and the eastern and western cluster sectors shown on Fig. [7] 
(in red and blue, respectively). Top: Gas surface brightness (left) and projected temperature profiles (right). Bottom left: ICM 
pressure profiles modelled to fit surface brightness and projected temperature. Bottom right: ICM pressure profiles re-scaled to 
the characteristic pressure of a self-similar cluster with mas s M-snn- The distribution of scaled pressure profiles measured in a 
representative sample of nearby galaxy clusters (REXCESS, lArnaud et al.ll2010l) has been overlaid in light blue. 



Tabl e 2. A2163 gas and total mass estimates performed from Yx mass proxy calibrated according to lArnaud et alJ {2010), 
IVikhlinin et alj d2009l) . hydrostatic equilibrium assumption in the overall cluster (Msoo.heX western (MsoonE.West) and eastern 
(A^500,HE,East) cluster sectors of Fig. [7] Confidence intervals in the Yx mass estimates include measurements uncertainties on gas 
mass and temperatures, and systematic uncertainties on the YX-M500 scaling relations. Uncertainties on hydrostatic mass estimates 
have been obtained from random realisations of parameters in the fitted density and temperature profile of Eq. (0 and (0. 
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6. Mass measurements 

The exceptionally high ICM temperature of A2163 and the ev- 
idence of an ongoing major merger event in its central region 
makes this cluster an interesting test case for cluster mass mea- 
surements. We first estimated the overall mass of A2163 from 
the Yx mass proxy, then investigated the validity of the hydro- 
static equilibrium assumption in the overall cluster and in spe- 
cific cluster sectors. 



6.1. Mass estimate from the Y x proxy 

Defined as the product of gas mass M g) 5oo and av erage tempera- 
ture kT , the Yx parameter has been proposed bv lKravtsov et al.l 
(2006) as a robust proxy for the total cluster mass. Related to 



the ICM thermal energy, the parameter has shown a tight corre- 
lation with cluster mass in mock X-ray observations of clusters 
evolving in a cosmological context, regardless of their dynami- 
cal state. 

We estimated the Yx parameter in A2163 iteratively, yield- 
ing the total cluster mass M500, radius ^00, gas mas M gi goo, 
and average temperature kT. We iterate about the Y x - M500 
scaling relation calibrated from hydrostatic mass estimates in 
a nearby cluster sample observed with XMM-Newton (A10), 
and compare with a similar relation obtained from Chandra data 
(IVikhlinin et al.ll2009l hereafter, V09). The gas masses M gj500 
were computed by integrating the ICM density profile used to 
model gas pressure as detailed in the previous section, while av- 
erage temperatures have been estimated within radii interval of 
[0. 15 - 0.75] r 50 o and [0.15 - 1] r 5 oo, following the Y x -M 50 o rela- 
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Fig. 9. A2163-A mass profiles computed assuming ICM hy- 
drostatic equilibrium. Profiles in the overall cluster region (in 
black), and the eastern and western cluster sectors illustrated in 
Fig G](in red and blue, respectively) are shown. The profiles have 
been derived from analytical ICM density and temperature dis- 
tributions fitted to the projected gas brightness and temperature 
profiles of Fig. [8] 



tions of A10 and V09, respectively. As shown in Table[2] the two 
Yx mass estimates obtained are consistent, with a lower estimate 
from XMM-Newton than from Chandra due to the difference in 
satellite calibrations. Our XMM-Newton cluster mass estimate 
from iteration about the V09 scaling relation is 7 per cent lower 
than that found by V09 using Chandra data. This mild differ- 
ence is likely to be related to the lower average temperature esti- 
mate of the present XMM-Newton analysis, possibly due resid- 
ual instrument cross-calibrations issues and small differences in 
Galactic absorption modelling. 

6.2. Mass estimates assuming ICM hydrostatic equilibrium 

Assuming hydrostatic equilibrium, the ICM density and temper- 
ature profiles used to model gas pressure have been used to de- 
rive the integrated cluster mass profile, shown as a grey shaded 
area in Fig. [9] As for the pressure measurements, this profile 
has been centred on the eastern X-ray peak of the Chandra map 
and estimated in the radius range 450-2100 kpc, excluding the 
central region and the crossing 'bullet' . As shown in Table [2] 
the total cluster mass estimate M500 obtained from this profile 
is inconsistent with the estimates obtained from the two XMM- 
Newton and Chandra Yx proxies (27 per cent and 19 per cent 
higher, respectively), possibly due to a departure of the ICM 
from hydrostatic equilibrium. 

To investigate the origins of this inconsistency, we derive 
two additional integrated cluster mass profiles using the density 
and temperature profiles from the eastern and western cluster 
sectors shown in Fig. [7] again using analytical models and as- 
suming ICM hydrostatic equilibrium. As shown in Fig. [9] the 
eastern and western cluster mass profiles significantly differ in 



shape. We observe that the western profile has a higher inte- 
grated mass than the eastern profile near the cluster centre but 
then flattens at larger radius, leading to a lower integrated mass 
estimate in the cluster outskirts. Similar to pressure profiles in 
Fig. [8] we further observe that the overall cluster mass profile 
appears to reflect the mass in the eastern cluster side. A2163 
being a complex cluster system, we cannot exclude the incon- 
sistency between eastern and western profiles to reflect a true 
anisotropy in the cluster mass distribution. The hypothesis of 
a strong mass excess beyond 1000 kpc in the eastern cluster 
side is however difficult to reconcile with the lack of any ob- 
servational evidence for prominent galaxy or da rk matter over- 
density in this cluster region -see paper I and iRadovich et alJ 
(120081) -. The inconsistency observed between eastern and west- 
ern profiles may instead show us that the shock heating of the 
cluster ICM induced by the accretion may not only have adi- 
abatically compressed the eastern cluster atmosphere, but also 
transiently taken gas away from hydrostatic equilibrium in this 
region. For this reason, the cluster mass estimates assuming hy- 
drostatic equilibrium in this region reported in Tableland hence 
in the overall cluster are likely to be overestimated. Assuming 
hydrostatic equilibrium in the western cluster side might instead 
provide with us an estimate of the cluster mass before the on- 
going merger event, Msoo,HE,West = 1-5 x 10 15 M , while the 
subcluster accretion might have raised the total cluster mass up 
to values provided by the Y x proxy: M 5 oo,y x - 1.9 x 10 15 M . 

7. Discussion and conclusions 

The A2163 system has been known to host a massive and elon- 
gated galaxy-cluster, A2163-A, and a smaller associated sub- 
cluster, A2163-B. We provided new mass measurements of the 
two companion subclusters, and estimate that with a mass of 
about M500 — 2.1 x 10 14 M , A2163-B is ten times less massive 
than its extended companion. As the ICM properties of A2163-B 
do not show any strong evidence of interaction with A2163-A, 
the two subclusters A2163-A are likely to be physically sepa- 
rated at a distance larger than main cluster radius (1500 - 2 Mpc). 
Since the two subclusters were suggested in Paper I to be con- 
nected by a galaxy bridge, they might form a bounded system 
observed along the line of sight, consistent with large separa- 
tion solutions in the two-body dynamical analysis detailed in 
Appendix. 

From evidence of a stripped cool core crossing the main clus- 
ter atmosphere separately from a nearby galaxy over-density, we 
infer that A2163 might have recently accreted a subcluster along 
its East- West elongation. We suggest that this merger event has 
shocked the main-cluster atmosphere, and adiabatically com- 
pressed the ICM behind the crossing cool core. This indication 
of shock heating and the evidence of a galaxy-gas separation 
lead us to infer that the subcluster has been accreted to a super- 
sonic velocity. Shock heating may further have transiently pulled 
cluster gas away from hydrostatic equilibrium, leading to large 
uncertainties on the cluster mass estimate from X-ray analysis. 
Assuming the Yx proxy to indicate the mass of the overall sys- 
tem A2163 and the western cluster sector to represent the ICM 
in hydrostatic equilibrium before subcluster accretion, we may 
infer a merger scenario between two subunits of mass-ratio 1 :4, 
leading to the present system with mass M500 - 1.9 x 10 M . 

A fast subcluster accretion with evidence of galaxy-gas sepa- 
rat ion was first seen i n IE 0657-56, the so-called 'bullet cluster' 
bv lMarkevitch et akl (120021) . It is interesting to note that both IE 
0657-56 and A2163 are exceptionally massive, favouring ram 
pressure stripping of the gas content of an incoming subclus- 
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ter due both to the high density of these accreting cluster atmo- 
spheres, and by exceptionally high collision velocities assuming 
free fall encounters. It is thus likely that the comparable mass of 
A2163 and IE 0657-56 has favoured similar merger scenarios. 
Moreover, both IE 0657-56 and A2163 exhibit similar global 
properties such as exceptionally high luminosity and gas tem- 
perature (about 15 and 12 keV for IE 0657-56 and A2163, re- 
spectively), and powerful emission from an extended radio halo. 
These common properties are likely point to the recent dissipa- 
tion of a large amount of kinetic energy through shock heating 
and turbulence. 

While presenting striking similarities related to a similar 
'bullet' accretion scheme, A2163 and IE 0657-56 show different 
morphologies that may point to some small differences between 
time scales and the initial conditions of the accretion process. 
The crossing 'bullet' is observed at a larger projected cluster ra- 
dius in the case of IE 0657-56 (500 kpc vs. 290 kpc as for IE 
0657-56 and A2163, respectively), and is thus better separated 
from the dense emission of the main cluster. Assuming both 'bul- 
let' motions to be roughly orthogonal to the line of sight, the pro- 
jected distances would represent physical distances to the cluster 
core, and the 'bullet' would be observed at an earlier point in its 
crossing of the main cluster in the case of A2163. Unlike for the 
IE 0657-56 'bullet', however, the A2163 'bullet' does not seem 
to be preceded by a clear shock front. This non detection might 
be the result of a less favourable merger geometry, with a larger 
angle separating the A2163 collision axis from the sky plane. As 
discussed in Paper I, the galaxy velocity distribution in the main 
cluster indeed shows a strong velocity gradient of ~ 1250 km/s 
along a NE/S W axis, which is much higher than the velocity off- 
set of ~ 600 km/s separating the two galaxy components of IE 
0657-56 (iBarrena et al.ll2002T) . 

As extensively studied in the textbook case of the 'bullet 
cluster', observing a high-velocity subcluster accretion within 
a massive cluster may provide interesting constraints on ICM 
physics, galaxy evolution, and dark matter properties. Better 
constraints on the dynamics of the exceptional merger event in 
A2163 would thus be obtained from deeper X-ray or SZ obser- 
vations, allowing us e.g. to constrain the 'bullet' velocities from 
precise pressure measurements at the cold front stagnation point, 
or to look for a shock front ahead of the 'bullet' direction of mo- 
tion. Moreover, the present observations might be compared to 
a high resolution weak lensing analyses, enabling us to verify 
if the bimodality of the galaxy distribution seen in Paper I re- 
flects the underlying dark matter morphology, and whether the 
galaxy-gas separation presently observed also corresponds to a 
segregation between the baryonic and dark matter content of the 
accreted subcluster. 
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Appendix: Two-body analysis of A2163A/A2163B 

We applied the two-body formalism dBeers et al.lll982l) to the 
system composed of A2163-A (main) and A2163-B (Northern 
subcluster). In this method, the two components are supposed 
to be coincident at t=0 and are moving apart or coming together 
for the first time in their history. In this framework, one can write 
the parametric solutions to the equations of motion, relating the 
elapsed time f, the separation R and the velocity V of the two 
components to the total mass of the system M, the development 
angle x, and the maximum expansion radius in the bound case 
(the asymptotic expansion velocity Vca, in the unbound case). In 
the case of bound radial orbits: 

R = ^.(l- C0&X ) (10) 
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Fig. 10. The sum of the virial masses of A2163-A and A2163-B as a function of the projection angle. The horizontal lines show 
the cluster mass estimate (full line) and errors (dashed lines) obtained in the present paper. The projected distance of 1 .25 Mpc and 
radial velocity difference of 200 km s _I are derived from Maurogordato et al. (2008). The two systems were at zero separation 10 
Gyr ago (left panel), 0.5 Gyr ago (middle panel), and 1.0 Gyr ago (right panel). The dotted curve corresponds to the Newtonian 
criterion for gravitational binding. Bound solutions are above this plot. 
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Table 3. Solutions of the two body model for the A2163-A- 
A2163-B system, for the three cases considered: before first ap- 
proach (10 Gyr), 0.5 Gyr and 1.0 Gyr after first approach. For 
each solution, the angle a between the line connecting the two 
components and the plane of the sky a, the spatial separation of 
the subclusters R, their separation at maximum expansion Rm 
and their relative velocity V are determined. 



3 \l/2 



\8GM/ 



(X ~ sin x) 



1/2 



V = 



sin x 



( 2GM \ 

\ R ln J (1 - COS X) 



In the unbound case: 



GM 

R = — (cosh x ~ 1) 
t = — (sinh^-^) 



V = V K 



sinh x 



(cosh x- 1) 



(11) 
(12) 

(13) 
(14) 
(15) 



The velocity and spatial separation are directly related to the ob- 
served line of sight velocity difference V, and projected separa- 
tion R p by: 



R = 



Rn 



-;V = 



sin a 



(16) 



where a is the angle between the axis joining the two units and 
the plane of the sky. When the two units are bound together, the 
system fulfils the Newton criterion: 



V r 2 R p < 2GM sin 2 a cos a 



(17) 



Specifying the angular separation, the velocity offset of the two 
systems, and the time, one can close the system and derive the 
evolution of the angle a as a function of these o bserved values 
and the development angle x (lOwers et al . 2009). This leads in 
the bound case to: 



tanff = 



V r t (1-cos^f) 2 



R p sinx(x - sinx) 
and similarly in the unbound case to: 

V r t (cOSh ^-l) 2 



tana- 



R p sinh^sinh^-^) 



(18) 



(19) 



These equations can be used in both cases to derive the total 
virial mass of the system as a function of the projection angle 
a with respect to the plane of the sky. Comparing to the current 
estimates o f the mass will allow us to assess which solutions are 
acceptable dBarrena et al.ll2002l) . 

We have taken R p = 1 .25 Mpc and V r = 200km/s (from 
Paper I), and set t to 10 Gyr (the age of the universe at z — 
0.2). We used the mass limits by summing mass estimates of 
the two components A2163-A and A2163-B as measured in the 
present paper and assuming M200 = F4M500. Fig.fTOlshows that 
in the pre-merger hypothesis, three bound solutions do exist: two 
bound incoming and one bound outgoing solutions (see Tablef3]l. 
The unbound solutions are excluded. 

We have also investigated the possibility that the first en- 
counter between both systems has already occurred, making sev- 
eral trials from t — 0.1 Gyr to lGyr. Recent encounters (f < 0.2 
Gyr) are not compatible with the system being bound. Taking 
t = 0.5 Gyr after the presumed encounter, the only possible so- 
lution is the bound outgoing solution. For higher values of f, one 
finds again three bound solutions (two incoming and one outgo- 
ing). 



